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We consider bilayer graphene in the presence of spin orbit coupling, to assess its behavior as a 
topological insulator. The first Chern number n for the energy bands of single and bilayer graphene 
is computed and compared, ft is shown that for a given valley and spin, n in a bilayer is doubled 
with respect to the monolayer. This implies that bilayer graphene will have twice as many edge 
states as single layer graphene, which we confirm with numerical calculations and analytically in 
the case of an armchair terminated surface. Bilayer graphene is a weak topological insulator, whose 
surface spectrum is susceptible to gap opening under spin-mixing perturbations. We also assess the 
stability of the associated topological bulk state of bilayer graphene under various perturbations. 
Finally, we consider an intermediate situation in which only one of the two layers has spin orbit 
coupling, and find that although individual valleys have non-trivial Chern numbers, the spectrum 
as a whole is not gapped, so that the system is not a topological insulator. 
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INTRODUCTION 

The study and characterization of the properties of 
topological insulators has sparked considerable interest 
Q] |2]. A topological insulator has a bulk energy gap 
separating the occupied electronic bands from the empty 
ones. However, the conduction and valence bands in the 
bulk are inverted with respect to their energy position 
in the vacuum. This necessarily results in the existence 
of localized surface states that cross the energy gap and 
are protected by time reversal symmetry. The electronic 
and magnetic properties of the surface states of a three- 
dimensional topological insulator are chiral, and are gov- 
erned by a two dimensional Dirac Hamiltonian. 

In a two dimensional topological insulator the chiral 
surface states become helical one dimensional edge states, 
and at a given edge, states with opposite spin orienta- 
tions transport charge in opposite spatial directions 3J. 
In the absence of time reversal symmetry breaking, the 
helicity of the edge states prevents backscattering be- 
tween states in the same edge and the Hall conductivity 
per spin becomes quantized, although with different sign 
for opposite spins. Note that this system does not vio- 
late time reversal symmetry because, when reversing the 
time direction, both the spin and the Hall conductivity 
per spin reverse. In this way the total Hall conductivity 
of the system vanishes. As pointed out by Thouless et al. 
[4] , the quantization of the Hall conductivity is a conse- 
quence of the topology of the band structure. Thouless 
et al. showed that the Hall conductivity should be an in- 
teger in units of e 2 /h, i.e., a topological invariant called 
the first Chern number. Because the Hall conductivity 
per spin is quantized, two-dimensional topological insu- 
lators are also called quantum spin Hall systems. The 
quantum spin Hall effect was theoretically predicted [5] 



to occur, and later experimentally observed [BJ, in HgTe 
quantum wells confined by CdTe barriers. 

Graphene, a two-dimensional carbon crystal [7], has 
also been proposed to be a two-dimensional quantum spin 
Hall system [SlUOj when the intrinsic spin orbit interac- 
tion is taken into account. Spin orbit coupling (SOC) in 
graphene opens a gap at the Dirac points and the system 
becomes a topological insulator. However, due to the n 
character of the graphene bands close to the Fermi en- 
ergy, the opened gap is proportional to the square of the 
intra-atomic spin orbit coupling constant divided by the 
energy difference between the s and p orbitals of graphene 
[Tflu"3"] . Since carbon is a very light atom, the spin orbit 
gap is very small, ~ 10 _3 meV. Consequently, the partic- 
ular transport properties of graphene as a topological in- 
sulator may only be observed in extremely clean samples 
and at extremely low temperatures. Graphene, neverthe- 
less, is a pedagogical toy model for analyzing properties 
of two-dimensional topological insulators [T] . 

In a Bernal stacked graphene bilayer, the conduction 
and the valence bands touch at two inequivalent points 
of the Brillouin zone and the system is a semimetal [14j . 
SOC in bilayer graphene also opens a gap in the band 
structure |15j . 

The spin-orbit interaction in bilayer graphene is larger, 
by about one order of magnitude, than in a single layer 
of graphene, due to the mixing of it and a bands by the 
interlayer hopping [TBJIT7]. Nevertheless, the spin-orbit 
induced gap in the bilayer is small and difficult to detect 
experimentally. As mentioned above, our main interest 
in this work is the use of graphene based structures to 
understand properties of topological insulators. In par- 
ticular, we use bilayer graphene as a model to analyze 
the coupling between two topological insulators. 

Our aim in this work is to study the topological prop- 



erties of a graphene bilayer in the presence of SOC. In 
graphene and bilayer graphene, the low energy proper- 
ties can be described by 2x2 Hamiltonians and for each 
momentum p, the wavefunctions have a spinor form. In 
undoped systems, the ground state of the system can 
be characterized by an unit vector field, h(p), that indi- 
cates, at each point of the reciprocal space, the expec- 
tation value of the orientation of the pseudospin. In the 
reciprocal space h has the form of a topological object. 
We find that in the case of graphene h has the form of 
a meron, and in the case of bilayer graphene it takes the 
form of a double vortex meron. Thus, we find that the 
first Chern number of bilayer graphene is twice that of 
the monolayer one, and consequently the number of edge 
states is also doubled. In bilayers, edge states are not 
Kramer-protected against backscattering, which give the 
topological insulating phase a weak character. We also 
analyze the stability of the insulating bulk topology in 
bilayer graphene with respect to bias voltage, staggered 
sublattice potential and trigonal warping effects. Finally, 
we study a bilayer graphene system in which only a sin- 
gle layer has SOC. We find that this system has a finite 
Chern number, but is a zero gap semiconductor for which 
no surface states are possible. 



GRAPHENE 

Carbon atoms in graphene are covalently bonded and 
arranged in a honeycomb lattice, which is composed of 
two triangular sublattices A and B. The low energy prop- 
erties in graphene are mainly determined by the 7r or- 
bitals. A tight-binding Hamiltonian with hopping 70 be- 
tween nearest-neighbors appropriately describes its band 
structure. In graphene the intrinsic SOC does not break 
the inversion symmetry of the honeycomb lattice and 
the electronic bands are spin degenerate. In addition, 
since the intrinsic SOC in graphene is a second order 
effect, the z-componcnt of the electron spin commutes 
with the Hamiltonian, and the bands can be indexed by 
s z = ±1, the up/down electron spin component perpen- 
dicular to the graphene layer. The corresponding tight 
binding Hamiltonian of graphene is 
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Here a^ Sz (bi tSz ) annihilates an electron on sublattice 
A(B) at site i and spin s z , t so is the next-nearest neigh- 
bor spin orbit hopping amplitude and Cj is either <ij or bj , 
depending wether the index j labels an A or B-sublatticc 
site, respectively. The factor z/y is +1 if the next-nearest 
neighbor hopping path rotates counterclockwise and -1 if 
it rotates clockwise. 

In undoped graphene, and for t so = 0, the conduction 
and valence bands touch at two inequivalent points of 



the Brillouin zone: K = ff (1,0) and K' = -f f (1,0), 
being a the lattice parameter of the triangular lattice. 
These are known as Dirac points. The main effect of the 
SOC in the electronic spectrum is the opening of a energy 
gap, 6\/3t so , at the Dirac points. Near these points the 
wave functions for each spin s z can be expressed via the 
k • P approximation |18j in terms of envelope functions 
^ Sz = [A a .(r),B..(r)] and ^_,„ = [<(r),B^(r)] for 
states near the K and K' points respectively. These 
wavefunctions satisfy the Dirac equations H Tz ^ Sz 'ip TzSz = 
eip Tz . Sz , where t z = ±1 specifies the Dirac points K and 
K' and 



H Tz ,a, = Vf(PxT z Vx +PyOy) + A so CJ z T z S 2 



(2) 



Here Vp= ^ Joa, A so = 3\/3t so and at are Pauli matri- 
ces representing the pseudospin degree of freedom corre- 
sponding to the two sites per unit cell of the graphene 
lattice. Note that p denotes the distance in momentum 
from the K and K' points. It is important to note that 
for a given s z , the Dirac mass terms in the Hamiltonian 
induced by SOC at K and K' have the same magnitude 
but different sign. 

Graphene with intrinsic SOC is a topological insulator 
with a finite spin Hall conductivity. For a given s z and 
valley index t z , the spin Hall conductivity as obtained 
from the Kubo formula has the form 
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where j(i) runs over occupied (empty) states, v v is the 
velocity operator in the v direction and the wavefunc- 
tions i\fu and energies £^ are obtained by diagonalizing 



Eq. (12|) for the corresponding t z and s z . In undoped 
graphene (Fermi energy crossing the Dirac points), the 
Hall conductivity for each valley takes the value 
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where h is the Planck constant. When summing over 
spins the total Hall conductivity of the system is zero, as 
it should be for a system with time reversal symmetry 
[TJ |2]. However, for each spin the Hall conductivity is 
quantized, although with opposite sign. This is the sig- 
nature of the quantum spin Hall effect. Note also that 
although for an isolated valley the Hall conductivity is a 
half integer in units of e 2 /h, the sum of the K and K' 
conductivities is quantized to integer multiples of e 2 /h, as 
it should be for a filled band of noninteracting electrons 
[HHH]. In an insulator, the value of the Hall conductivity 
in units of e 2 jh is related to the first Chern number [4] 
of its bandstructure. The Chern number corresponding 
to a 2x2 Hamiltonian 
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is related to the number of times the unit sphere is cov- 
ered by the unit vector h(p) when p runs over the whole 
reciprocal space [U El [20] . This number takes the form 



4tt 



d 2 p(d Px h x d Py h) ■ h, 



(6) 



where h(p) depends on r z and s z and ±e(p) are the energy 
eigenvalues, n is a topological invariant, the Pontyagin 
index of the mapping h(p) [20;. Physically, the vector 
h(p) represents the expectation value of the orientation 
of the pseudospin associated with the wavefunctions of 
the Hamiltonian, Eq. (pi). 

The vector h defines the topology of the band structure 
and may be written in the form 
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Here the valley 



with h z (p) = s z T z A so /yJv 2 F p 2 + A 2 so 
index t z = ± defines the right and left vorticity of the 
topological structure and 9 is the azimuthal angle made 
by the momentum vector p. At asymptotically large mo- 
mentum, h z vanishes, while in the cortex core we have 
h z = t z s z . This implies that there are four flavors 
of topological objects which are usually referred to as 
merons [2"Tll2"2"] . since they are essentially half skyrmions. 
The four possible merons are illustrated in Fig. [TJ The 
Chern number index corresponding to the field h takes 
the form 
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Thus, the topological charge corresponding to each Dirac 
point is ±1/2 depending on the sign of the electron spin. 
When summing over the contribution from both Dirac 
points the Chern number in graphene is ±1. This result 
follows from the fact that, topologically, a meron has half 
the winding number of a skyrmion, that is, the topolog- 
ical object formed by the two Dirac points together. 



Helical Edge States 

The most spectacular consequence of the existence 
of an insulator with a topologically non trivial band 
structure is the appearance of gapless conducting states 
at interfaces where the topology of the band structure 
changes. In two dimensional quantum spin Hall sys- 
tems, helical edge states appear at the surfaces of the 
material. At each edge there exists a pair of one di- 
mensional channels with opposite spins that propagate 
in opposite directions. In graphene these states have 
been obtained numerically by diagonalizing nanoribbons 
terminated with different geometries, zigzag [TUl |23j or 
armchair [23]. Here we describe analytically how helical 




FIG. 1: (Color online) Four flavors of merons, see Ref. 1211 
These are vortices that are right or left handed depending on 
the valley index t z . The topological charge, ±1/2, is deter- 
mined by the sign of the spin s z . 



one dimensional channels appear at an armchair termi- 
nated edge. 

We consider an edge running along the y-direction and 
the vacuum region is defined by the x < condition. 
The graphene armchair termination consists of a line of 
A-B dimers, so it is natural to have the wave function 
amplitude vanish on both sublattices at X = 0. To do 
this we must admix valleys [2"5] and require 



A Sz (x = 0) 4- A' Sz (x = 0) = and 
B Sz {x = 0) + B' s (x = 0) = 0. 



(9) 



Note that the boundary conditions should be satisfied by 
each spin separately. We consider solutions with momen- 
tum p y along the surface and energy E inside the gap, 
\E\ < A so , where p x = ihn. The general solution has the 
form 
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with Hvfk — \ I &-1o — E 2 + Vpp 2 and 
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The boundary conditions, Eq. (19]), imply that tan£ = 

tan£'. This, together with the definition of k, gives the 

solutions 



E = s z vfP v and hvpn = A s 



(12) 



Since the graphene flake exists in the x > region, nor- 
malizability of the wavefunction implies that n, should be 



denoted by e^' (p), are 
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FIG. 2: (Color online) Tight-binding band structure of an 
armchair terminated wide nanoribbon. Shadow regions rep- 
resent bulk states. 



positive for any spin. Therefore, quasiparticles with pos- 
itive spin, s z > 0, move in the positive y direction and 
have and energy E — vpp y . On the contrary, quasipar- 
ticles with s z < move in the negative y direction and 
have energy E = —vpp y . On the opposite surface, for 
which the wavefunction exists for x < 0, normalizability 
implies that p x = —ihn. In this case E = —s z vfP v , and 
the edge states with positive (negative) velocity in the 
y-direction have s z < (s z > 0). In Fig. [2] we plot the 
band structure of a wide armchair terminated nanorib- 
bon with a SOC t so = O.OI70. 



GRAPHENE BILAYERS 
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These energies are independent of s z and r z [26 . The 

(2) 
eigenvalues e± describe two strong interlayer coupling 

(2) (2) 

bands with energies e + > 71 and e_ < —71. These 
bands do not touch at the Dirac points, and correspond 
to wavefunctions mostly localized at sites A\ and B2, 
which form strong dimers [57]. The eigenvalues e± de- 
scribe low energy bands. Performing the usual low energy 
approximation vpp -C 71 , the dispersion of these bands 

where m 
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can be approximated by £± ; = ±yj £^ + Af , 
is an effective mass induced by the interlayer hopping, 
171 = 7i/2i>fi. The corresponding low energy eigenstates 
are mostly localized at sites B\ and Ai- All of these 
states are degenerate in the spin s z and valley t z indices. 
For a given s z and r z , the Hall conductivity can be ob- 
tained numerically by plugging the eigenvalues and eigen- 
functions of the Hamiltonian ( 13 1 into the expression pi. 
For charge-neutral bilayer graphene, only negative energy 
bands arc filled, and the Hall conductivity takes the form 
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When summing over the valleys, the Hall conductivity 
per spin is twice that of graphene. Interestingly, this 
is the same result one would expect for the case of two 
decoupled layers, although the eigenfunctions are com- 
pletely different from those of the coupled bilayer. As in 
the case of a monolayer, time reversal symmetry dictates 
opposite a^y Tz upon reversing the spin, so that the total 
Hall conductivity of the system is zero. 



Bilayer graphene consists of two stacked graphene 
sheets. In bilayer graphene there are four sites per unit 
cell, which we label A\, B\ and A2 and B2 in the the 
first and second layer respectively. We consider the so- 
called Bernal stacking, commonly found experimentally, 
in which site Bi is exactly on top of the sublattice A\ . 
Interlayer coupling is modeled by hopping amplitude 71 
between these two sites in each unit cell. The low en- 
ergy properties of this model are well described for each 
spin s z by the following Hamiltonian acting on the four- 
component spinor (Ai, B\, A 2 , B 2 ) [T4"j . 
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Band structure topology 

The value of the Hall conductivity can be understood 
from the topology of the band structure through the first 

Chern number. Because of the gap ~ 271 between the 

(2) 
high energy dimer bands £ ± , these cannot change their 

topology by adding a weak SOC. Therefore, only the 
gapless £j_ bands acquire non-trivial topology that con- 
tributes to a finite Chern number n. To compute n, it is 
then sufficient to consider the effective 2x2 low energy 
Hamiltonian that describes the low energy bands of the 
system in the limit e <C 71 [13] . 



where H TztSs _ is the monolayer Hamiltonian of Eq. (12]), ai 
and Ti are the Pauli matrices for the sublattice and layer 
degree of freedom respectively and 7o is the unit matrix 
in the layer subspace. The four energy bands of H BG , 
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FIG. 3: (Color online) Four flavors of double vortex merons. 
These are vortices that are right or left handed depending on 
the valley index r z . The topological charge ±1 is determined 
by the sign of the spin. 



FIG. 4: (Color online) Band structure of a wide armchair 
terminated bilayer graphene ribbon for a single spin. The 
(green) shaded areas represent the continuum bulk states. 



This Hamiltonian acts on the two-component spinor 
(Bi, A 2 ). This Hamiltonian describes the low energy ef- 
fective coupling between carbon atoms in different layers 
which are not directly connected by tunneling. Their 
coupling arises as a result of virtual transitions through 
the high energy dimer states that have been integrated 
out. From this Hamiltonian and the relation of Eq. ^ 
we get the following expression for the unit vector field 
h(p): 



h = 
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As in the case of 



A? 



with h z (p) = A so r z s z /y 

a monolayer, the valley index t z determines the vorticity 
of the field h(p) and the angle 8 is the azimuthal angle of 
momentum p. At large momentum p, the z-component 
of the field h vanishes and the field is confined to be in 
the x — y plane. As in the case of the monolayer, there 
are four flavors of these topological objects, shown in Fig. 
pi The Chern number of the corresponding fields h(p) is 



n = 2- 



Att 
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p dp 



Thus, the Chern number corresponding to each valley is 
±1 depending on the sign of the electron spin. It is in- 
teresting to note that, although these objects have twice 
the Pontryagin charge of merons (which are often thought 
as half a skyrmions), they are not skyrmions. Far away 
from the center of a skyrmion, i.e., at p — > oo, the orien- 
tation of the vector field h(p) is constant. This is not the 
case of bilayer graphene (see Fig. [3]), although, as in a 
skyrmion, the Chern number is unity for each t z ,s z . We 
call these objects double vortex merons. When summing 



over valleys the Chern number in bilayer graphene is ±2 
(depending on the spin), in agreement with the results 
obtained for the Hall conductivity. 



Edge states 

According to the bulk-boundary correspondence rule 
PQ, and since the Chern number of bilayer graphene is 
twice that of the monolayer, the number of helical sur- 
face states should also be double. In order to analyze 
the edge states, we diagonalize the tight-binding Hamil- 
tonian corresponding to a wide armchair terminated bi- 
layer graphene ribbon. In Fig. HI we plot its band struc- 
ture as function of the wavevector k y along the ribbon for 
a given electron spin s z . As in the case of a monolayer 
of graphene, the Hamiltonian commutes with s z and the 
band structure is degenerate in spin. The bulk bands 
are indicated by shaded regions. We identify two pairs 
of parallel edge states inside the gap. The states with 
positive velocity are located on the opposite ribbon edge 
from the states with negative velocity. This can be seen 
by analyzing the eigenfunctions of the low-energy Hamil- 
tonian, Eq. (16 1. As a result of the helicity of the edge 



states, their location changes to the opposite edge when 
the spin of the carriers changes sign. 

Bilayer graphene edge states can be understood as 
bonding/antibonding combinations of the surface states 
of the two constituent graphene layers. In the absence of 
tunneling between the graphene layers, the surface states 



have the form of Eq. ( 10 1 with £, £', E and K determined 
by Eq. (Ill and Eq. (12). For finite 71, surface states 



from layers 1 and 2 become coupled. To leading order 
in 71, the energies and wavefunctions of the graphene 




FIG. 5: (Color online) Schematic diagram of the edge states 
in a bilayer graphene system. At the edge of the sample there 
are four channels: two with spin up and moving in one direc- 
tion, and two with spin down an moving in the opposite direc- 
tion. Dotted and dashed lines represent channels formed by 
bonding/antibonding combinations of individual layers edge 
channels. 



bilayer Hamiltonian Eq. (13 1, defined for x < 0, are 
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where ± correspond to bonding/antibonding states. 
These expressions are valid when the SOC A s0 is much 
larger than 71 /4. In other situations, bulk states should 
reduce the bonding/antibonding energy separation. 

By inspection of Eq. ( p0| , we see that at each edge 
of the system there exist four one dimensional channels: 
two with spin up moving in one direction and the other 
two with spin down and moving in the opposite direc- 
tion, see Fig. [5j The two terminal electrical conductance 
in this model is 4e 2 /h. The injection of charge current 
through edge states would result in antisymmetric spin 
accumulation at the edges. In a four terminal geometry 
the spin conductances should be quantized. 

However, in bilayer graphene, backscattering between 
channels at the same edge moving in opposite directions 
with opposite spins is not forbidden, in contrast to the 
case of monolayer graphene. Indeed, Kramers's theorem 
ensures degeneracy of time-reversed pairs in the absence 
of time reversal symmetry breaking, preventing backscat- 
tering between energy degenerate edge states with oppo- 
site spins at high symmetry points of the Brillouin zone, 
such as the T-point (for which k = —k modulo a recip- 
rocal lattice vector) [Tj. While in monolayer graphene 
edge states branches cross at the T-point k y = 0, bilayer 
graphene edge state branches cross away from k y = 0, see 



Fig. H As a result, SOC perturbations will in general 
split their degeneracy, and open a gap in the edge state 
dispersion around zero energy. Therefore, this model is 
not, strictly speaking, a topological insulator, since there 
exist single-body effects that are symmetric under time 
reversal that can induce backscattering between edge 
channels, and make the edges insulating. This happens 
for any two-dimensional insulator with an even number 
N K of edge state pairs, related by time reversal symme- 
try (k y ,s z — > —k y ,—s z ), at any given energy inside the 
gap (Nk = 2 for the case at hand). An odd number of 
pairs Nk, like in monolayers or trilayers with SOC, im- 
plies that there exists at least one of those pairs crossing 
at a high symmetry point, and are therefore protected 
from scattering by Kramer's theorem. A Z 2 topological 
invariant of the form 



v = iV^mod 2, 



(21) 



is defined pQ to differentiate strict topological insulators 
(dubbed 'strong'), v = 1, from 'weak' topological insu- 
lators, with v = and Nk 7^ 0. Within this classifi- 
cation, bilayer graphene is a weak topological insulator. 
The 'weakness' is in relation to the class of time-reversal- 
invariant perturbations that can open a gap in the edge 
states, in this case general spin-mixing perturbations. 

Stability of the topological insulating phase 

The non-trivial insulating phase created by the instrin- 
sic SOC may be destroyed by Rashba SOC or other per- 
turbations if they are strong enough. This happens not 
by inducing backscattering between edge states, but by 
band reconnection, which changes the topology of the 
bulk bandstructure back to that of a conventional insu- 
lator jS]. It is a generic possibility both in weak and 
strong topological insulators. 

We now consider the destruction of bilayer graphene's 
topological insulator phase by a symmetric Rashba SOC, 
a staggered sublattice mass term and a voltage bias be- 
tween the two layers. We also discuss the effect of the 
trigonal warping on the topological insulating phase. 
1) The Rashba spin orbit coupling term arises due to an 
electric field perpendicular to the bilayer plane or from 
the interaction with a substrate, 



H R = \r% ® (t z o- x 



(22) 



Here s x and s y are the Pauli matrices in the spin sub- 
space and Xr is the intensity of the coupling. This term 
violates mirror symmetry about the planes. In the pres- 
ence of this term the z-component of the electron spin 
is not conserved, but there is a region of values of the 
coupling, < Xr < A so , for which the ground state 
is adiabatically connected [HI HH] to the topological insu- 
lating phase. For values of the coupling A^ > A so the 



energy gap closes and the electronic structure is that of 

a zero gap semiconductor with quadratically dispersive 

bands. 

2) A staggered sublattice potential of the form 



H r 



A„7o ® cr z 



(23) 



stabilizes the band insulator phase that competes with 
the topological insulating phase. This term is typically 
zero in graphene but would be present for a similar boron 
nitride film. In the effective 2x2 Hamiltonian this term 
appears as a mass term independent of the spin and val- 
ley, X v cr z . In the presence of thi s term, the form of the 
unit vector field h(p) of Eq. (17) is unaffected, but now 
the z-component takes the form 



h z (p) = 



t-^so7~zS Z i ^v 



(A so T z s z + X v ) 2 + ^y 
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FIG. 6: (Color online ) D ensity plot in the reciprocal space of 
the integrand of Eq. ( 18 1 in the absence (a) and the presence 
(b) of trigonal warping. The center of each square corresponds 
to the original Dirac point at K. Cold colors represent regions 
where the integrand is bigger. The parameters used in these 
plots are 71=0.170, A so = 10 _J 7o and in (b) 73=0.0570. 



The Chern number associated with this vector field h(p) 

is 



n = -t z sign(A so T z s z + X v ) 



(25) 



For values A„ < A so the Chern number is —s z (equal for 
both valleys) and the system is a topological insulator. 
For values X v > A so we have n = —t z and the total Chern 
number per spin is zero when summing over valleys. In 
this case the bilayer graphene becomes a band insulator. 
3) A bias voltage between the layers gives a contribution 



to the bilayer Hamiltonian, Eq. (13), of the form 



Hi 



bias 



= VT Z 



<7 . 



(26) 



This term acts equally on both sublattices and on both 
spin orientations. In the effective 2x2 low energy Hamil- 



tonian Eq. (16), this term has exactly the same form as 



the staggered sublattice potential Va z . Thus, the condi- 
tion for the existence of a topological insulating phase in 
the bilayer is V < A so . 

In the last two cases, when X v > A so or V > A so , the 
system becomes a conventional insulator with vanishing 
Hall and spin Hall conductances. The two inequivalent 
valleys, considered independently, have nonzero topolog- 
ical charge with opposite sign, n — — t z , and a valley de- 
pendent Hall effect may occur (2S]. The nonzero Chern 
number of each valley induces chiral modes in topological 
domain walls in graphene [29] and bilayer graphene |30j . 
4) Trigonal warping. The inclusion in the Hamiltonian 
of bilayer graphene of a weak direct A2-B1 hopping term 
73 <C 71 introduces in the 2x2 effective Hamiltonian a 
term of the form 



H„ 



V3(T z p x O- x - PyO-y), 



(27) 



with 113 = ^ajs/h. This term produces a trigonal warp- 
ing in the band structure, which stretches the isoenergy 



lines along the directions = 0, |7r and |7r for the val- 
ley K and along the directions cf> — g7r, it and |7r for 
the valley K'[37]. For A so =0, the trigonal warping pro- 
duces a dramatic change in the band structure. For a 
given valley, instead of two parabolic bands touching at 
the Dirac point, there are now four Dirac points, one at 
the center (K or K'-point), and three others that occurs 
at finite momentum and in the directions <p mentioned 
above. The SOC opens gaps at the four Dirac points, 
making the system insulating. To analyze the nature 
of the insulating phase we compute the Chern number. 
The warping term breaks the s ymm etric form of the unit 
vector field h(p) given in Eq. ([17]) and prevents an ana- 
lytical calculation of the Chern number. By integrating 
numerically Eq. (13J) we obtain that the Chern number 
is not affected by the trigonal warping. In the absence 
of trigonal warping the main contribution to the integral 
comes from an annulus around each Dirac point where 
dh z /dp is maximum, see Fig. [61a) (dark region). When 
the trigonal warping effects are taken into account, the 
main contribution to the Chern number comes from the 
regions around the three new Dirac points, Fig. p[ 



SPIN-ORBIT PROXIMITY EFFECT 

In this Section we analyze the case of a graphene bi- 
layer in which only one of the layers has a finite intrinsic 
SOC A so , whereas the other layer has zero coupling. This 
could happen in a bilayer in contact with a strong SOC 
metal [31], or in a bilayer where heavy atoms or molecules 
have been deposited on top of one of the two layers [3T]. 

When only one of the layers is affected by the SOC, 
the 2x2 effective Hamiltonian takes the form 



Ht z ,s 



<-^so7~ z S z 
1m 






FIG. 7: (Color online) Low energy band structure of a 
graphene bilayer where only one of the two layers has an finite 
SOC A so . Left and right panels correspond to the bands near 
the K and K' points, respectively. Bands for states with spin 
up s z = 1 (down s z = — 1) are represented with solid (dashed) 
lines. The Fermi energy is indicated by an dotted horizontal 
line. The SOC opens a gap A so with a different sign for each 
valley and each spin orientation. 
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7~zPxPy 



<j£28) 



and the unit vector field h has almost the same form as 
in Eq. (17). Thus, the Chern number of this system is 



equal to that of a graphene bilayer with intrinsic SOC in 
both layers, 

n = -s z . (29) 

Although from the value of the Chern number the system 
appears to be a topological insulator, this is not the case: 
it turns out it is not a true insulator at all. From the 
diagonalization of the Hamiltonian of Eq. (28) we get 
the band dispersion 



J 7~zS ; 



e± 



±\ A 



im 2 



(30) 



This band structure, shown in Fig. [7J corresponds to 
an indirect zero band gap semiconductor. The absence 
of SOC in one of the layers breaks the inversion sym- 
metry, making the bands non-degenerate under indepen- 
dent inversion of valley and spin. When both valleys 
and/or both spins are considered, there is no gap in the 
spectrum. However, for each spin and valley there is an 
energy gap A so at p = 0, and the system is thus an in- 
sulator from an optical point of view. Moreover, if the 
system is free from perturbations that mix the two val- 
leys and the two spins, it will also behave as an effec- 
tive electronic insulator. This has an important conse- 
quence on the existence of surface states in the bilayer 
structure. In zigzag edge terminations for which it is 
not necessary to admix valleys to satisfy the boundary 



FIG. 8: (Color online) Tight binding band structure of a 
zigzag (a) and a armchair (b) terminated bilayer graphene 
nanoribbon in which there is SOC only in one of the layers. 
These results correspond to s z = l. Results for s z =-l are ob- 
tained by changing k to —k. In the numerics 71 = O.I70 and 
A so = 3\/3 x 10~ 2 7o. The width of the armchair and zigzag 
nanoribbons is 151a and 120\/3ffl respectively. In the arm- 
chair case all the states are extended across the ribbon. In 
the zigzag one the states in the center of the Brillouin zone 
are localized at the edges. 



conditions [25|. states coming from different valleys can 
be treated independently and dispersive surface localized 
channels appear in the spectrum, see Fig. ^[a). However, 
for armchair terminated ribbons where, in order to satisfy 
the boundary conditions, the wavefunctions have to be a 
linear combination of both valleys, no localized surface 
states appear in the band structure, see Fig. [8[b). The 
absence of a full energy gap and the lack of surface states 
in some boundary terminations are clear indicators that 
a graphene bilayer with SOC in only one layer is not a 
quantum spin Hall system. 



CONCLUSION 

In this work we reviewed the description of single layer 
graphene with SOC as a topological insulator in terms 
of the first Chern number, which arises naturally in the 
computation of the Hall conductivity. We showed then 
that the Chern number per spin in bilayer graphene is 
two, twice that of the monolayer. This doubling is re- 
flected also in the number of topological surface states, 



as a consequence of the bulk-surface correspondence rule. 
By numerically computing the spectrum of finite size 
samples, we checked that bilayer systems have twice as 
many edge states as a monolayer. This was furthermore 
confirmed analytically in the case of surfaces with arm- 
chair termination. The fact that the total Chern number 
per spin is even means that the bilayer system is a weak 
topological insulator, susceptible to gap opening if the 
system is subject to spin-mixing perturbations (such as 
a weak Rashba SOC). We also assessed the general sta- 
bility of the bulk topological insulating state of bilayer 
graphene with respect to Rashba SOC, a staggered sub- 
lattice potential, interlayer bias, and trigonal warping. 
The first three perturbations compete with the intrinsic 
SOC and, when sufficiently large, spoil the inverted gap 
property that is crucial to making the system a topolog- 
ical insulator. Finally, we examined a bilayer graphene 
system in which only one layer has intrinsic SOC. Al- 
though in this system individual valleys have non-trivial 
Chern numbers, the spectrum as a whole is not gapped, 
so that the system is not a topological insulator. 

Note added: While this manuscript was in the final 
stage of preparation, a manuscript by Cortijo, Grushim 
and Vozmediano appeared|33j which studies the Chern 
Simons coefficients in graphene and bilayer graphene us- 
ing an effective action formalism. 
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